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Abstract 

This paper presents the results of studies on the outlet boundary conditions for turbulent internal flow 
simulations. Several outlet boundary conditions have been investigated by applying the National 
Combustion Code (NCC) to the configuration of a LM6000 single injector flame tube. First of all, very 
large eddy simulations (VLES) have been performed using the partially resolved numerical simulation 
(PRNS) approach, in which both the nonlinear and linear dynamic subscale models were employed. 
Secondly, unsteady Reynolds averaged Navier-Stokes (URANS) simulations have also been performed 
for the same configuration to investigate the effects of different outlet boundary conditions in the context 
of URANS. Thirdly, the possible role of the initial condition is inspected by using three different initial 
flow fields for both the PRNS/VLES simulation and the URANS simulation. The same grid is used for all 
the simulations and the number of mesh element is about 0.5 million. 

The main purpose of this study is to examine the long-time behavior of the solution as determined by 
the imposed outlet boundary conditions. For a particular simulation to be considered as successful under 
the given initial and boundary conditions, the solution must be sustainable in a physically meaningful 
manner over a sufficiently long period of time. 

The commonly used outlet boundary condition for steady Reynolds averaged Navier-Stokes (RANS) 
simulation is a fixed pressure at the outlet with all the other dependent variables being extrapolated from 
the interior. The results of the present study suggest that this is also workable for the URANS simulation 
of the LM6000 injector flame tube. However, it does not work for the PRNS/VLES simulation due to the 
unphysical reflections of the pressure disturbances at the outlet boundary. This undesirable situation can 
be practically alleviated by applying a simple unsteady convection equation for the pressure disturbances 
at the outlet boundary. The numerical results presented in this paper suggest that this unsteady convection 
of pressure disturbances at the outlet works very well for all the unsteady simulations (both PRNS/VLES 
and URANS) of the LM6000 single injector flame tube. 

1.0 Introduction 

Large eddy simulation (LES) and very large eddy simulation (VLES) of reacting turbulent internal 
flows are critically important for the accurate modeling of the mixing and combustion processes occurring 
in the combustors. Recently, we have developed an approach, called the partially resolved numerical 
simulation (PRNS), which aims at bringing out the dynamically important unsteady large and very large 
scale turbulent structures in the numerical simulation, but only using computer resources typically 
required by the unsteady Reynolds averaged Navier-Stoke approach (URANS). As a result, the 
PRNS/VLES approach is quite practical for engineering application. Its physical fidelity is higher than 
that of the RANS, while its demand on the computing resources is significantly lower than the demand by 
the LES. 
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In our previous efforts (Refs. 1 to 4), either the periodic boundary condition or the extrapolation 
boundary condition was imposed at the outlet of the computational domain; mainly for the purpose of 
assessing the fundamentals and/or demonstrating the short-time solution of the PRNS/VLES approach. 
However, it is well known that the outlet boundary condition is critically important to the long-time 
solution of numerical simulation of turbulent flows (Refs. 5 to 6). The goal of the present effort is to 
identify a workable but simple outlet boundary condition (BC) for the long-time PRNS/VLES solution of 
turbulent internal flows. Our criteria for a workable boundary condition are set as follows: It must lead to 
a physically meaningful numerical solution; and this meaningful solution must be sustainable over a 
sufficiently long period of time. 

In the current assessment effort, we have chosen the convective type of boundary condition (Ref. 7) 
as our baseline. Other types of outlet boundary conditions such as the extrapolation boundary condition 
and the fixed pressure boundary condition are also used for comparison. The numerical simulations were 
performed using the configuration of a LM6000 single injector flame tube, as its characteristic flow 
features are representative of the flow features typically occurring in the practical combustors. The same 
grid of about 0.5 million elements was used in all of the calculations. 

Both PRNS/VLES and URANS numerical simulations were carried out. Two subscale models 
(nonlinear and linear) of PRNS/VLES were applied. Two Reynolds stress models (nonlinear and linear) 
of URANS were also used. For each type of simulation (PRNS/VLES or URANS), different outlet 
boundary conditions (i.e., the convective BC, the extrapolated BC, and the fixed pressure BC) were 
applied for comparison. Furthermore, three different initial conditions (i.e., the nonlinear RANS solution, 
the linear RANS solution and the static flow field) were employed for each type of simulation to examine 
the effect of the initial conditions. 

The numerical results are presented in terms of the time history of flow variables at four locations 
along the centerline of the flame tube; the instantaneous distribution of flow variables at a center plane, 
and the variation of flow variables along the centerline. The effects of the initial condition and the 
turbulence model are examined in the context of the imposed outlet boundary condition for either the 
PRNS/VLES or the URANS simulation. 

2.0 Unsteady Outlet Boundary Conditions 

A brief description of the unsteady outlet boundary condition and its implementation in the NCC will 
be given in this section. 

2.1 Unsteady Convective Boundary Condition (BC) 

The following unsteady convective boundary condition (Ref. 7) will be applied at the outlet 
boundary: 


dt 


+ U 


fljj H 

dn 


Jb 


= 0 , 


$£ = W, V, IV,/*, p, h, k, 8 


( 1 ) 


Where (/ represents a dependent flow variable, for example, the velocity components u, v, w, the mass 
fraction of species f k , the gauge pressure p and the specific enthalpy /z, the subscale turbulent kinetic 
energy and its dissipation rate k and s . The outlet boundary is indicated by the subscript b and its unit 
outward normal is denoted by n. U is a global “convective” velocity out of the boundary surface, and its 
magnitude is determined, at any instant, by the requirement of the global mass conservation of the entire 
computational domain. 
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2.2 Implementation of Unsteady Convective BC in the NCC 

Consistent with the solution algorithm for the interior flow field (Ref. 8), Eq. (1) is also solved in a 
time-accurate manner via the same dual time step scheme in which the convergence of the inner (or the 
pseudo time) loop is achieved through the application of a 4 stage Runge-Kutta scheme, i.e., 


§\b = 

$b = 
$b = 
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The residual is defined as 


R 




‘ + ^ $b 


\+u l 


Wtb 

8n 


(3) 


Here, n represents the pseudo time. Within each pseudo time, m (= 0, 1,2, 3) is related to the stage 
number of the Runge-Kutta scheme, and p denotes the real time. Ax is the pseudo time step and At is the 
real time step. At every real time p , Eq. (2) is executed for pseudo time n ( 1 — » go ) until the residual R is 
reduced to a prescribed order of magnitude. This procedure at the boundary is synchronized with the 
solver for the interior field to provide the necessary boundary information at each and every Runge-Kutta 
stage. 


2.3 Extrapolation Boundary Condition 

The extrapolation BC has the following form: 


f 

V 


8n 


Jb 


= 0 


(4) 


3.0 PRNS/VLES of LM6000 Single Injector Flame Tube 

LM6000 is a General Electric low NOx emission gas turbine combustor. We have performed several 
types of numerical simulations for a single injector flame tube to evaluate the effects of several outlet 
boundary conditions. A highly swirling jet (a lean methane-air mixture) is injected from a circular inlet 
into a rectangular duct. The inlet pressure and temperature are about 6 atmospheres and 644 K. The inlet 
flow variables are specified by using the mean profiles from the experiment. The Reynolds number based 
on the inlet axial velocity and the inlet jet diameter is about 3,200,000. The following two figures depict 
the computational domain and the numerical grids on two perpendicular center planes. The total number 
of grid points is about 495,000, it is noted here that this same grid is used in this study for all simulations 
(PRNS/VLES, URANS and RANS). 
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LM6000, 495,000 elements 
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In this section, we will present the results from very large eddy simulation using the PRNS approach. 
The resolution control parameter (RCP) was set at a value of about 0.333. Two subscale models 
(nonlinear and linear) were applied. Three outlet boundary conditions (i.e., the unsteady convective BC, 
the extrapolation BC and the fixed pressure BC) were tested. Furthermore, for a given outlet boundary 
condition, three different initial conditions (i.e., the nonlinear RANS solution, the linear RANS solution 
and the static flow field) were used in its assessment. 


3.1 Nonlinear Subscale Model 

PRNS with the nonlinear subscale model is fundamentally different from the traditional LES 
approach. Here, the interactions between the resolved large scale turbulence and the unresolved small 
scale turbulence are accounted for not just by the eddy viscosity, but also, explicitly in the filtered 
transport equations, by the turbulence source terms originated from the nonlinear part of the subscale 
model. So far, these additional turbulence source terms have not been considered in the existing LES 
approaches (see Ref. 1). 

3.1.1 Results Using Unsteady Convective BC for Pressure Only 

In this case, the unsteady convective outlet BC is only applied to the gauge pressure. The rest of the 
dependent flow variables at the outlet are determined by extrapolating from the interior point (i.e. 
extrapolation BC). This is consistent with the observation that, for subsonic viscous internal flow 
simulation, the information on pressure at the outlet boundary is always needed for properly maintaining 
the global mass conservation. 

The results are presented in three parts: the time history of velocity components and gauge pressure at 
four locations along the centerline; the instantaneous contour plots of flow variables at a center plane; and 
the instantaneous centerline flow variable profiles. 

3. 1.1.1 Time History 

The time history of velocity components and gauge pressure were recorded at four centerline 
locations: x = 0.015, 0.05, 0.10 and 0.2. From which we may examine the development of turbulent 
fluctuations and perform further spectra analysis. The results are presented with respect to three different 
initial conditions: the nonlinear RANS solution, the linear RANS solution, and the static flow field. 
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Initial condition: Linear RANS solution 
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Time histry at Probe 1 Time histry at Probe 2 
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Initial condition: Static flow field 

Time histry at Probe 1 Time histry at Probe 2 
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3. 1.1.2 Contour at Center Plane 


The snapshots of instantaneous flow field at a center plane are presented for the time step number 
60,000, which is about 100 flow-through time (defined as the ratio of the length of the combustor to the 
inlet centerline axial velocity). From which we may examine the flow structures. The results are presented 
with respect to three different initial conditions: the nonlinear RANS solution, the linear RANS solution, 
and the static flow field. 


Initial condition: Nonlinear RANS solution 
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Initial condition: Linear RANS solution 
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PRNS, non-lineae, RCP=0.3 
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Initial condition: Static flow field 
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Extrapolated BC at exit for others 




turb_k 


n 


6.92E+01 
643E+01 
5.93E+01 
544E+01 
4 94E+01 
445E+01 

3 95E+01 
346E+01 
2.97E+01 
247E+01 
1.98E+01 
1 48E+01 
9.89E+00 

4 94E+00 


J 


Y 


PRNS, non-lineae, RCP=0.3 
2nd= 0.0, 4th=0.01 , CFL = 1 , dt = 4e-6 
0.5 M elements, 70,000 time steps 
Starting from fresh run 
Refined convective BC at exit forpg 
Extrapolated BC at exit for others 



mu 

1.06E-02 

9.85E-03 

9.09E-03 

8.33E-03 

7.58E-03 

6.82E-03 

6.06E-03 

5.30E-03 

4.55E-03 

3.79E-03 

3.03E-03 

2.27E-03 

1.52E-03 

7.58E-04 



PRNS, non-lineae, RCP=0.3 
2nd= 0.0, 4th=0.01 , CFL = 1 , dt = 4e-6 
0.5 M elements, 70,000 time steps 
Starting from fresh run 
Refined convective BC at exit forpg 
Extrapolated BC at exit for others 





P9 

3.35E+03 

1.33E+03 

-6.95E+02 

-2.71E+03 


■ -4.74E+03 

■ -6.76E+03 

■ -8.78E+03 
-1.08E+04 
-1.28E+04 

| -148E+04 

■ -1.69E+04 
-1.89E+04 

■ -2.09E+04 
™ -2.29E+04 



Y Vorti city Magnitude 

3.00E+04 
2.77E+04 
2.54E+04 
2.31 E+04 
2.08E+04 
1 .85E+04 
1.62E+04 
1.39E+04 
1.16E+04 
9.30E+03 
7.00E+03 
4.70E+03 
2.40E+03 
1.00E+02 



PRNS, non-lineae, RCP=0.3 
2nd= 0.0, 4th =0.01 , CFL = 1 , dt = 4e-6 
0.5 M elements, 70,000 time steps 
Starting from fresh run 
Refined convective BC at exit forpg 
Extrapolated BC at exit for others 
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3. 1.1.3 Centerline Variables for Three Different Initial Conditions 

The instantaneous centerline flow variables (axial velocity u , Mach number, subscale turbulent 
kinetic energy k , effective viscosity \i T , gauge pressure and vorticity magnitude) at the time step 60,000 
are presented here with three different initial conditions (the nonlinear RANS solution, the linear RANS 
solution, and the static flow field). From which we may examine the flow variation along the centerline. 
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3.1.2 Results Using Extrapolation BC 

In this case, all dependent flow variables at the outlet, including the gauge pressure, are extrapolated 
from the interior point. This is reminiscent of the so called perfectly non-reflecting boundary condition 
which might be used without significant numerical problem for short-time simulations (Ref 6). 

The results are presented in three parts: the time history of velocity components and gauge pressure at 
four locations along the centerline; the instantaneous contour plots of flow variables at a center plane; and 
the instantaneous centerline flow variable profiles. 

3. 1.2.1 Time History 

The time history of velocity components and gauge pressure were recorded at four centerline 
locations: x = 0.015, 0.05, 0.10 and 0.2. From which we may examine the development of turbulent 
fluctuations and perform further spectra analysis. The simulations have been carried out using three 
different initial conditions: the nonlinear RANS solution, the linear RANS solution, and the static flow 
field. They all lead to the conclusion that, in the long run, the numerical solution can not be sustained in a 
physically meaningful manner. In the following, only the results using the linear RANS solution as the 
initial condition are presented to illustrate the development of the calculation. 

The time histories at all four locations show that the flow variables are not approaching their 
statistically stationary mean values even after a long-time simulation. 
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200 


Time histry at Probe 3 
PRNS,RCP=0.3, Non-linear 
2nd=0.0, 4th=0 .01 , dt=4.0e-6 
Starting from RANS linear solution 
Extrapolated BC at Exit 
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3. 1.2.2 Contour at Center Plane 


The snapshots of instantaneous contour of flow variables at a center plane are presented for two time 
instants, namely, the 20,000 and the 55,000 time steps. It is evident that the flow field has undergone a 
dramatic change over this period of time. We consider that the flow field at time step 20,000 is still 
reasonable, but, as time goes by, the dynamically important flow structures are not sustained, in fact, the 
solutions are becoming physically unreasonable. 



PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01, cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 20,00 time step 
Starting from RANS_L solution 
Extrapolated BC 




u 

1.00E+02 

9.12E+01 

8.19E+01 

7.27E+01 

6.34E+01 

5.41E+01 

4.49E+01 

3.56E+01 

2.63E+01 

1.71E+01 

7.80E+00 


-1.47E+00 


-1.07E+01 

-2.00E+01 




PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01, cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 55,000 time step 
Starting from RANS_L solution 
Extrapolated BC 




u 

1 .90E+02 
1 .74E+02 
1 .58E+02 
1.42E+02 
1 .25E+02 
1 .09E+02 
9.31E+01 
7.69E+01 
6.08E+01 
4.46E+01 
2.85E+01 
1 .23E+01 


-3.85E+00 


-2.00E+01 




PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01, cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 20,00 time step 
Starting from RANS_L solution 
Extrapolated BC 


u 


Mach 
2.80E-01 
2.60E-01 
2. 4 IE-01 
2.21E-01 
2.01E-01 
1.81E-01 
1.61E-01 
1.41E-01 
1.21E-01 
1.01E-01 

8.1 IE-02 

6.1 IE-02 
4.12E-02 
2.13E-02 




PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01, cfl = 1 .0, dt = 4.E-06 
0.5 M elements, at 55,000 time step 
Starting from RANS_L solution 
Extrapolated BC 




Mach 
3.93E-01 
3.65E-01 
3.37E-01 
3.10E-01 
2.82E-01 
2.55E-01 
2.27E-01 
1.99E-01 
1.72E-01 
1 .44E-01 
1.16E-01 
8.88E-02 
6.12E-02 
3.35E-02 
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Vorticity Magnitude 
3.00E+04 
2.77E+04 
2.54E+04 
2.31E+04 
2.08E+04 
1.85E+04 
1.62E+04 
1.39E+04 
1.16E+04 
9.30E+03 
7.00E+03 
4.70E+03 
2.40E+03 
1.00E+02 



Y Vorticity Magnitude 

n 3.00E+04 
2.77E+04 
2.54E+04 
2.31E+04 
2.08E+04 


PRNS, RCP=0.3, non-linear model ui 1.85E+04 

2nd= 0.0, 4th=0.01 , cfl = 1 .0, dt = 4.E-06 1.62E+04 

0.5 M elements, at 55,000 time step 1.39E+04 

Starting from RANS_L solution 1.I6E+04 

Extrapolated BC 9.30E+03 


7.00E+03 



P 

PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01, cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 20,00 time step 
Starting from RANS_L solution 
Extrapolated BC 



PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01, cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 20,00 time step 
Starting from RANS_L solution 
Extrapolated BC 




turb_k 
3.56E+00 
3.31E+00 
3.05E+00 
2.80E+00 
2.54E+00 
2.29E+00 
2.04E+00 
1 .78E+00 
1.53E+00 
1.27E+00 
1 .02E+00 
7.63E-01 
5.09E-01 
2.54E-01 




PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01, cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 55,000 time step 
Starting from RANS_L solution 
Extrapolated BC 


turb_k 
— 7.50E-01 
■ 6.96E-01 
B 6.43E-01 
I 5.89E-01 
5.36E-01 
4.82E-01 
4.29E-01 
3.75E-01 
3.21E-01 
2.68E-01 
2.14E-01 
1.61E-01 
1.07E-01 
5.36E-02 




PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01, cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 20,00 time step 
Starting from RANS_L solution 
Extrapolated BC 


mu 

I 7.77E-04 
7.24E-04 
6.7 IE-04 
6.18E-04 
5.64E-04 
5.1 IE-04 
4.58E-04 
4.04E-04 
3. 5 IE-04 
2.98E-04 
2.45E-04 
1.9 IE-04 
1 .38E-04 
8.49E-05 




PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01, cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 55,000 time step 
Starting from RANS_L solution 
Extrapolated BC 


mu 

m 4.52E-05 
B 4.42E-05 
4.32E-05 
4. 2 IE-05 
4.1 IE-05 
4.0 IE-05 
IJ 3. 9 IE-05 
P 3. 8 IE-05 
3. 7 IE-05 
3. 6 IE-05 
3. 5 IE-05 

■ 3.40E-05 
3.30E-05 
3.20E-05 
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Y 


j=L 


PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01, cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 20,00 time step 
Starting from RANS_L solution 
Extrapolated BC 


pg 

I 1.60E+05 
1.58E+05 
1 .55E+05 
1.53E+05 
1.50E+05 
1.48E+05 
1.45E+05 
1 .43E+05 
1.41E+05 
1.38E+05 
1.36E+05 

■ 1.33E+05 
1.31E+05 
1.29E+05 



Y 



PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01 , cfl = 1 .0, dt = 4.E-06 
0.5 M elements, at 55,000 time step 
Starting from RANS_L solution 
Extrapolated BC 




pg 

-1.54E+05 

-1.57E+05 

-1.59E+05 

-1.61E+05 

-1.64E+05 

-1.66E+05 

-1.69E+05 

-1.71E+05 

-1.74E+05 

-1.76E+05 

-1.78E+05 

-1.81E+05 

-1.83E+05 

-1.86E+05 



3.1.3 Results Using Fixed Pressure at the Outlet 

In this case, fixed gauge pressure is imposed at the outlet while the rest of the variables are 
extrapolated from the interior point. The nonlinear RANS solution is used as the initial condition. 

The calculation crashes after 7650 time step. Here we present the instantaneous contour of flow 
variables at the center planes for two time instants: 2550 time step and 7650 time step. 


3. 1.3.1 Contour at Center Planes 

Snapshots of flow variable contours are shown at tow perpendicular center planes. Figures below on 
the left side are the snapshots at the time step 2,550, and figures on the right side are the snapshots taken 
at the time step 7650, i.e., when the simulation is about to crash. 

The results indicate that, although the fluctuating turbulent flow field can be established over a 
relatively short time period, eventually, it can not survive the impact of pressure disturbances reflected 
from the outlet boundary having a fixed pressure. Very large inflow appears in the outlet region and the 
corresponding Mach number can exceed 1.4. Obviously, the flow structures are totally unphysical right 
before the crash of the calculation. 



PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01 , cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 2550 time step 
starting from steady RANS solution 
Using fixed pg at exit 


u 

8.91E+01 
7.84E+01 
6.77E+01 
5.70E+01 
4.63E+01 
3.56E+01 
249E+01 
1 42E+01 
3.50E+00 
-7.20E+00 
-1.79E+01 
-2.86E+01 


-3 93E+01 
-5 00E+01 




PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01, cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 7650 time step 
starting from steady RANS solution 
Using fixed pg at exit 


u 

891E+01 
7.84E+01 
6.77E+01 
5 70E+01 
4 63E+01 
3.56E+01 
2.49E+01 
1 42E+01 
3.50E+00 
-7.20E+00 
-1.79E+01 
-2 86E+01 


-3.93E+01 

-5.00E+01 
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z 


z 



PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01 , cfl = 1 .0, dt = 4.E-06 
0.5 M elements, at 2550 time step 
starting from steady RANS solution 
Using fixed pg at exit 


n 8 91E+01 
7 84E+01 
6.77E+01 
5.70E+01 
4.63E+01 
3.56E+01 
2.49E+01 
1.42E+01 
"1 3.50E+00 
-7.20E+00 
-1.79E+01 


y 


-2.86E+01 
-3.93E+01 
-5 00E+01 




PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01 , cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 7650 time step 
starting from steady RANS solution 
Using fixed pg at exit 




8.91E+01 
7.84E+01 
6.77E+01 
5.70E+01 
4.63E+01 
3.56E+01 
2.49E+01 
1 42E+01 
3.50E+00 


-7.20E+00 


-1.79E+01 


■ -2.86E+01 
® -5.00E+01 




PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01 , cfl = 1 .0, dt = 4.E-06 
0.5 M elements, at 2550 time step 
starting from steady RANS solution 
Using fixed pg at exit 


Mach 
M 3.13E-01 
I 2.90E-01 
■ 2.68E-01 
*\ 2.46E-01 
2.24E-01 
U 2.02E-01 
1.79E-01 
1.57E-01 
1.35E-01 
1.13E-01 
9.04E-02 
6.82E-02 
4.60E-02 
2.38E-02 




PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01, cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 7650 time step 
starting from steady RANS solution 
Using fixed pg at exit 


Mach 


u 


1 39E+00 
1.30E+00 
1 20E+00 
1.10E+00 
9.97E-01 
8.97E-01 
798E-01 
698E-01 
5.99E-01 
4.99E-01 
4 00E-01 
3.00E-01 
2.00E-01 
1.01E-01 




PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01 , cfl = 1 .0, dt = 4.E-06 
0.5 M elements, at 2550 time step 
starting from steady RANS solution 
Using fixed pg at exit 


Mach 
3.13E-01 
2.90E-01 
2.68E-01 
2.46E-01 
2.24E-01 
2.02E-01 
I 1.79E-01 
■ 1.57E-01 
1.35E-01 
1.13E-01 
9.04E-02 

■ 6.82E-02 
4.60E-02 
2.38E-02 




PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01 , cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 7650 time step 
starting from steady RANS solution 
Using fixed pg at exit 


Mach 




1.39E+00 
1 30E+00 
1.20E+00 
1.10E+00 
9.97E-01 
8.97E-01 
7.98E-01 
6.98E-01 
5.99E-01 
4.99E-01 
4.00E-01 
3.00E-01 
2.00E-01 
1.01E-01 
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PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01 , cfl = 1 .0, dt = 4.E-06 
0.5 M elements, at 2550 time step 
starting from steady RANS solution 
Using fixed pg at exit 


Vorticity Magnitude 


3.00E+04 

2.77E+04 

2.54E+04 

2.31E+04 

2.08E+04 

1.85E+04 

1.62E+04 

1.39E+04 

1.16E+04 

9.30E+03 

7.00E+03 

4.70E+03 

2.40E+03 

1.00E+02 



Y 

a 


PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01 , cfl = 1 .0, dt = 4.E-06 
0.5 M elements, at 7650 time step 
starting from steady RANS solution 
Using fixed pg at exit 


Vorticity Magnitude 

n 3.00E+04 
2.77E+04 
2.54E+04 
A 2.31E+04 
2.08E+04 
1.85E+04 
1.62E+04 
1.39E+04 
1.16E+04 
9.30E+03 
7.00E+03 
4.70E+03 
J 2.40E+03 
1.00E+02 




PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01 , cfl = 1 .0, dt = 4.E-06 
0.5 M elements, at 2550 time step 
starting from steady RANS solution 
Using fixed pg at exit 




pg 

3.93E+03 
1 .92E+03 
-8.41E+01 
-2.09E+03 
-4.10E+03 
-6.1 0E+03 
-8.11E+03 
-1.01E+04 
-1.21E+04 
-1.41E+04 
-1.61E+04 
-1.81E+04 
-2.0 IE +04 
-2.22E+04 




PRNS, RCP=0.3, non-linear model 
2nd= 0.0, 4th=0.01, cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 7650 time step 
starting from steady RANS solution 
Using fixed pg at exit 


I 


pg 

1.80E+05 
1.58E+05 
1.37E+05 
1.15E+05 
9.38E+04 
7.23E+04 
5.08E+04 
2.92E+04 
7.69E+03 
-1 38E+04 
-3.54E+04 


u 


-5.69E+04 

-7.85E+04 

-1.00E+05 



3.2 Linear Subscale Model 

By now, it is clear that PRNS with the nonlinear subscale model and using the unsteady convective 
outlet boundary condition for the gauge pressure can successfully simulate the flow in a LM6000 single 
injector fame tube starting from different initial conditions. In the following, we turn our attention to 
PRNS with the linear subscale model. 

A linear subscale model is a pure eddy viscosity model that is used in all the existing LES type of 
simulation, in which the effects of the unresolved small scale turbulence on the resolved large scale 
turbulence are solely accounted for via the eddy viscosity. In our previous assessment effort focusing on 
the fully developed turbulent pipe flows (Ref. 1), it has been demonstrated that the linear subscale model 
is not adequate for simulations of low Reynolds number turbulent pipe flows, because the turbulent 
fluctuations can not be sustained in the simulation over a long period of time. 

3.2.1 Results Using Unsteady Convective BC for Pressure Only 

Here, we present the results using the linear subscale model for the simulation of a LM6000 single 
injector flame tube. Both of the initial and the outlet boundary condition are the same as that applied in a 
case of using the nonlinear subscale model. That is, we apply the unsteady convective BC for the gauge 
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pressure while the rest of the dependent variables are extrapolated from the interior; and the initial 
condition is the nonlinear RANS solution. 

The results indicate that, although the short-time solution looks reasonable, this simulation does not 
provide the long-time solution, more specifically, the calculation crashes after 18,000 time steps. 


3.2.1. 1 Contour at Center Plane 


Snapshots of flow variable contours are shown at a center plane for two instances: the time step 
10,000 (figures on the left below) and the time step 18,000 (figures on the right below). At time step 
10,000 (about 20 flow-through time), the flow structures appear to be reasonable, but the subscale 
turbulent kinetic energy k and eddy viscosity p r are too small, they are about two orders-of-magnitude 
smaller than their counterparts when using the nonlinear subscale model. At time step 18,000 (about 36 
flow-through time), the flow structures near the outlet become unphysical, namely, large amount of flow 
rushing into the domain, with its corresponding Mach number approaching one. The calculation crashes 
soon after. We attribute the failure of the liner subscale model in this case to the fact that, by its nature, it 
can not account for the rotational and anisotropy effects. Hence, it is not adequate to use the linear 
subscale model for simulating the high swirling flow occurring in this LM6000 single injector flame tube. 


PRNS, linear model, RCP=0.3, 

2nd= 0.0, 4th=0.01, cfl=1, dt=4e-6 
0.5 M elements, at 10,000 time step 
Starting from RANS_NL solution 
refined convective BC at exit for pg 
Others extrapolated 



u 

n 1.06E+02 
9.66E+01 
8.69E+01 
7.72E+01 
■ 6.75E+01 
I 5.78E+01 
I 4.80E+01 
I 3.83E+01 
2.86E+01 
1.89E+01 
9.16E+00 
■ -5.60E-01 
-1.03E+01 


-2.00E+01 



PRNS, linear model, RCP=0.3, 

2nd= 0.0, 4th=0.01, cfl=1, dt=4e-6 
0.5 M elements, at 18,000 time step 
Starting from RANS_NL solution 
refined convective BC at exit for pg 
Others extrapolated 





u 

1.06E+02 

9.66E+01 

8.69E+01 

7.72E+01 

6.75E+01 

5.78E+01 

4.80E+01 

3.83E+01 

2.86E+01 

1.89E+01 

9.16E+00 


-5.60E-01 


-1.03E+01 


-2.00E+01 



PRNS, linear model, RCP=0.3, 

2nd= 0.0, 4th=0.01 , cfl=1 , dt=4e-6 
0.5 M elements, at 10,000 time step 
Starting from RANS_NL solution 
refined convective BC at exit for pg 
Others extrapolated 




Mach 
3.13E-01 
2.91 E-01 
2.68E-01 
2.46E-01 
2.24E-01 
2.02E-01 
1.79E-01 
1.57E-01 
1.35E-01 
1.13E-01 
9.05E-02 
6.83E-02 
4.60E-02 
2.38E-02 



PRNS, linear model, RCP=0.3, 

2nd= 0.0, 4th=0.01 , cfM , dt=4e-6 
0.5 M elements, at 18,000 time step 
Starting from RANS_NL solution 
refined convective BC at exit for pg 
Others extrapolated 



n 


Mach 
3.33E-01 
3.09E-01 
2.85E-01 
2.62E-01 
2.38E-01 
2.14E-01 
1.91 E-01 
1.67E-01 
1.43E-01 
1.20E-01 
9.58E-02 
7.21 E-02 
4.84E-02 
2.48E-02 
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Mach contour at exit 
PRNS, linear model, RCP=0.3 
2nd= 0.0, 4th=0.01, cfl=1, dt=4e-6 
0.5 M elements, at 10,000 time step 
Starting from RANS_L solution 
refined convective BC at exit for pg 
Others extrapolated 




Mach 

n l .0000E-01 
9.2308E-02 
8.4615E-02 
7.6923E-02 
6.9231 E-02 
6.1538E-02 
5.3846E-02 
4.61 54E-02 
3.8462E-02 
3.0769E-02 
2.3077E-02 

■ 1.5385E-02 
7.6923E-03 
0.0000E+00 


Mach contour at exit 
PRNS, linear model, RCP=0.3 
2nd= 0.0, 4th=0.01, cfl=1, dt=4e-6 
0.5 M elements, at 18,000 time step 
Starting from RANS_L solution 
refined convective BC at exit for pg 
Others extrapolated 




Y 


Mach 


I 


9.1460E-01 
8.4424E-01 
7.7389E-01 
7.0354E-01 
6.3318E-01 
5.6283E-01 
4.9248E-01 
4.2212E-01 
3.51 77E-01 
2.8141 E-01 
2.1 1 06E-01 
1.4071 E-01 
7.0354E-02 
0.0000E+00 


PRNS, linear model, RCP=0.3, 

2nd= 0.0, 4th=0.01, cfl=1, dt=4e-6 
0.5 M elements, at 10,000 time step 
Starting from RANS_NL solution 
refined convective BC at exit for pg 
Others extrapolated 


tL 




Vorticity Magnitude 
3.00E+04 
2.77E+04 
2.54E+04 
2.31 E+04 
2.08E+04 
1.85E+04 
1.62E+04 
1.39E+04 
1.16E+04 
9.30E+03 
7.00E+03 
4.70E+03 
2.40E+03 
1.00E+02 



PRNS, linear model, RCP=0.3, 

2nd= 0.0, 4th=0.01, cfl=1 , dt=4e-6 
0.5 M elements, at 18,000 time step 
Starting from RANS_NL solution 
refined convective BC at exit for pg 
Others extrapolated 


Vorticity Magnitude 
3.00E+04 
2.77E+04 
2.54E+04 
2.31 E+04 
2.08E+04 
1.85E+04 
1.62E+04 
1.39E+04 
1.16E+04 
9.30E+03 
7.00E+03 
4.70E+03 
2.40E+03 
1.00E+02 



turb k 


PRNS, linear model, RCP=0.3, 

2nd= 0.0, 4th=0.01 , cfl=1 , dt=4e-6 
0.5 M elements, at 10,000 time step 
Starting from RANS_NL solution 
refined convective BC at exit for pg 
Others extrapolated 


CL 



PRNS, linear model, RCP=0.3, 

2nd= 0.0, 4th=0.01, cfl=1, dt=4e-6 
0.5 M elements, at 18,000 time step 
Starting from RANS_NL solution 
refined convective BC at exit for pg 
Others extrapolated 





turb_k 

5.00E-01 
4.62E-01 
4.23E-01 
3. 85 E-01 
3.46E-01 
3.08E-01 
2. 69 E-01 
2.31 E-01 
1 .92 E-01 
1 .54 E-01 
1.15E-01 
7.69E-02 
3.85E-02 
0.00E+00 
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PRNS, linear model, RCP=0.3, 

2nd= 0.0, 4th=0.01, cfl=1, dt=4e-6 
0.5 M elements, at 10,000 time step 
Starting from RANS_NL solution 
refined convective BC at exit for pg 
Others extrapolated 


I 


mu 

8.00E-05 
7.58E-05 
7.15E-05 
6.73E-05 
6.31 E-05 
5.88E-05 
5.46E-05 
5.04E-05 
4.62E-05 
4.19E-05 
3.77E-05 
3.35E-05 
2.92E-05 
2.50E-05 



PRNS, linear model, RCP=0.3, 

2nd= 0.0, 4th=0.01, cfl=1 , dt=4e-6 
0.5 M elements, at 18,000 time step 
Starting from RANS_NL solution 
refined convective BC at exit for pg 
Others extrapolated 


Y 



mu 

8.00E-05 
7.58E-05 
7.15E-05 
6.73E-05 
6.31 E-05 
5.88E-05 
5.46E-05 
5.04E-05 
4.62E-05 
4.19E-05 
3.77E-05 

3.35E-05 

2.92E-05 

2.50E-05 





pg 

5.37E+03 

3.20E+03 

1.02E+03 

-1.15E+03 

-3.33E+03 

-5.50E+03 

-7.68E+03 

-9.86E+03 

-1.20E+04 

-1.42E+04 

-1.64E+04 

-1.86E+04 

-2.07E+04 

-2.29E+04 



PRNS, linear model, RCP=0.3, 

2nd= 0.0, 4th=0.01, cf 1= 1 , dt=4e-6 
0.5 M elements, at 10,000 time step 
Starting from RANS_NL solution 
refined convective BC at exit for pg 
Others extrapolated 


n 


PRNS, linear model, RCP=0.3, 

2nd= 0.0, 4th=0.01, cfl=1, dt=4e-6 
0.5 M elements, at 18,000 time step 
Starting from RANS_NL solution 
refined convective BC at exit for pg 
Others extrapolated 



pg 


n 


1.23E+04 
9.78E+03 
7.26E+03 
4.74E+03 
2.23E+03 
-2.93E+02 
-2.81 E+03 
-5.33E+03 


-7.85E+03 


-1.04E+04 


-1.29E+04 

-1.54E+04 

-1.79E+04 


-2.04E+04 



3.3 Concluding Remarks 

The unsteady pressure convective BC appears to be very suitable for PRNS simulations with 
nonlinear subscale model. It is also robust with respect to very different initial conditions. It is noticed 
that the solutions with different initial conditions are different from each other near the inlet region where 
no or less turbulent fluctuations are observed due to the specified “laminar” inflow condition; however, 
the solutions become statistically equivalent in regions away from the inlet, where the turbulence 
becomes fully developed. The solution does sustain itself in a physically meaningful manner over a long 
period of time. 

Both extrapolation BC and fixed pressure at the outlet are not suitable for large or very large eddy 
simulations. Even though a “reasonable” solution can be established, they can only last over a short 
period of time. In the long run, the solution will become unphysical. 

In addition, the numerical studies indicate that the linear subscale k-s model does not work very 
well for high swirling flows. The subscale eddy viscosity is vanishing towards the value of the laminar 
viscosity, eventually, the simulation can not survive for long in the high Reynolds number environment. 
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4.0 URANS of LM6000 Single Injector Flame Tube 

In this section, we present the results of unsteady RANS simulations. Two Reynolds stress models 
(nonlinear and linear) have been applied. Three outlet boundary conditions (i.e., the unsteady convective 
BC, the extrapolation BC and the fixed pressure BC) have been assessed. Furthermore, three different 
initial conditions (i.e., the nonlinear RANS solution, the linear RANS solution and the static flow field) 
have been used 

4.1 Nonlinear Model 

URANS with the nonlinear model is fundamentally different from the standard k-e model (Ref. 9). 
Here, the interactions between the turbulent mean flow and the entire spectrum of the turbulent 
fluctuations are accounted for not only by the eddy viscosity, but also by the explicit turbulence source 
terms in the filtered transport equations. These additional source terms are due to the nonlinear part of the 
model, and they do not originate from the standard k - 8 model. 

4.1.1 Results Using Unsteady Convective BC for Pressure Only 

Similar to the PRNS simulations in Section 3.1.1, the unsteady convective outlet BC is only applied 
to the gauge pressure. All the other dependent flow variables at the outlet are provided by extrapolating 
from the interior point. 

The results are presented in three parts: the time history of velocity components and gauge pressure at 
four locations along the centerline; the instantaneous contour plots of flow variables at a center plane; and 
the instantaneous centerline flow variable profiles. 

4.1. 1.1 Time History 

The time history of velocity components and gauge pressure are recorded at four centerline locations: 
x = 0.015, 0.05, 0.10 and 0.2. From which we may examine the temporal development of the filtered 
variables. The results are presented with respect to three different initial conditions: the nonlinear RANS 
solution, the linear RANS solution, and the static flow field. 

Initial condition: Nonlinear RANS solution 


Time history at Probe 1 
80 r- URANS, Non-linear model 
- 2nd=-0.01 , 4th=0.05, dt=4.0e-6 

_ Starting from RANS_NL solution 

Refined convective BC for pg, others extrapolated 

60 - _ 


£ 40 - 
> 


Time history at Probe 2 
60 i- URANS, Non-linear model 

- 2nd=-0.01, 4th=0.05, dt=4.0e-6 

- Starting from RANS_NL solution 

- Refined convective BC for pg, others extrapolated 
40 - 


> 


20 - 


u 


v 

w 


20 - 


0 - 


-20 - 
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5000 

10000 

15000 

20000 

-™ 0 

5000 

10000 

15000 

20000 



Time step 





Time step 
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M ‘A 


Time history at Probe 3 
20 i- URANS, Non-linear model 

- 2nd=-0.01, 4th=0.05, dt=4.0e-6 

- Starting from RANS_NL solution 

- Refined convective BC for pg, others extrapolated 

10 - u 

V 

w 


0 - 


Time history at Probe 4 
50 i- URANS, Non-linear model 
I 2nd=-0.01 , 4th=0.05, dt=4.0^6 
: Starting from RANS_NL solution 

40 - Refined convective BC for pg, others extrapolated 


30 - 



> 


3 io - 


u 

V 


w 


-10 - 


-20 I 1 ■ I ■ ■ ■ ■ I I I 1:11,1 

0 5000 10000 15000 20000 

Time step 


o - 


-io - 


-20 r | 1 
0 


I I I I I I I I I I I I I I I I 

5000 10000 15000 20000 

Time step 


Time history at Probe 1 



Time history at Probe 2 



Time history at Probe 3 Time history at Probe 4 
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Initial condition: Linear RANS solution 


Timehistry at Probe 1 
180 r URANS, Non-linear model 

2nd=-0.01, 4th=0.05, dt=4.0e-6 
1 60 I- Starting from RANS linear solution 


I Refined convective BC for pg, others extrapolated 


140 - 
120 - 
J00 - 

- 80 - 

3 

60 - 
40 - 
20 - 

0 - 


5000 


10000 

Time step 


15000 


20000 


60 r 


40 - 


-20 - 


Timehistry at Probe 2 

URANS, Non-linear model 

2nd=-0.01, 4th=0.05, dt=4.0e-6 

Starting from RANS linear solution 

Refined convective BCfor pg, others extrapolated 


_ (ki Lb 


10000 

Time step 


20 r 


£ 

> 0 


-10 - 


Timehistry at Probe 3 

URANS, Non-linear model 

2nd=-0.01, 4th=0.05, dt=4.0e-6 

Starting from RANS linear solution 

Refined convective BCforpg, others extrapolated 






5000 


10000 

Time step 


15000 


20000 


50 r 


40 - 


30 - 


20 - 


-20 


Time histry at Probe 4 

URANS, Non-linear model 

2nd=-0.01, 4th=0.05, dt=4.0e-6 

Starting from RANS linear solution 

Refined convective BC for pg, others extrapolated 




0 - 


10000 

Time step 
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v, w u,v, w 


Initial condition: Static flow field 


Timehistry at Probe 1 

URANS, Non-linear model 

2nd=-0.01, 4th=0.05, dt=4.0e-6 

Starting from static flow field 

Refined convective BCfor pg, others extrapolated 


60 - 


40 - 




J i i i i I i i i i I i i i i I i i i i I 

0 5000 10000 15000 20000 

Time step 


Timehistry at Probe 2 



Timehistry at Probe 3 



Time step 


Timehistry at Probe 4 
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4.1. 1.2 Contour at Center Plane 


The snapshots of instantaneous contour of flow variables at a center plane are presented for the time 
step 20,000, which is about 40 through-flow time. From which we may examine the flow structures. The 
results are presented with respect to three different initial conditions: the nonlinear RANS solution, the 
linear RANS solution, and the static flow field. 


Initial condition: Nonlinear RANS solution 


Y 


URANS, non-lineae model 
2nd= -0.01 , 4th=0.05, CFL = 1 , dt = 4e-6 
0.5 M elements, 20,000 time steps 
Starting from steady RANS_NL solution 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 




u 

n 1.08E+02 
9.78E+01 
8.80E+01 
7.82E+01 
6.84E+01 
5.86E+01 
4.87E+01 
| 3.89E+01 
2.91E+01 
1.93E+01 
9.46E+00 

■ -3.61E-01 
-1.02E+01 
-2.00E+01 


Y 


URANS, non-lineae model 
2nd= -0.01 , 4th=0.05, CFL = 1 , dt = 4^6 
0.5 M elements, 20,000 time steps 
Starting from steady RANS_NL solution 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 



n 


Mach 
3.1 0E-01 
2.88E-01 
2.66E-01 
2.44E-01 
2.22E-01 
2.00E-01 
1.77E-01 
1.55E-01 
1.33E-01 
1.1 IE-01 
8.88E-02 
6.66E-02 
4.45E-02 
2.24E-02 



Y 


URANS, non-lineae model 
2nd= -0.01 , 4th=0.05, CFL = 1 , dt = 4^6 
0.5 M elements, 20,000 time steps 
Starting from steady RANS_NL solution 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 



l 


turb_k 
1.95E402 
1 .81E+02 
1 .67E+02 
1 .53E+02 
1 .39E+02 
1 .26E+02 
1.12E+02 
9.76E+01 
8.37E+01 
6.97E+01 
5.58E401 
4.18E401 
2.79E401 
1 .39E+01 



Y 


URANS, non-lineae model 
2nd= -0.01 , 4th=0.05, CFL = 1 , dt = 4e-6 
0.5 M elements, 20,000 time steps 
Starting from steady RANS_NL solution 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 



I 


mu 

7.23E-02 

6.72E-02 

6.20E-02 

5.68E-02 

5.17E-02 

4.65E-02 

4.13E-02 

3.62E-02 

3.10E-02 

2.58E-02 

2.07E-02 

1.55E-02 

1.04E-02 

5.19E-03 



Y 


URANS, non-lineae model — . 

2nd= -0.01 , 4th=0.05, CFL = 1 , dt = 4e-6 x 

0.5 M elements, 20,000 time steps 
Starting from steady RANS_NL solution 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 


D 


pg 

8.55E+03 

6.37E+03 

4.19E+03 

2.02E+03 

-1.60E+02 

-2.34E+03 

-4.51E+03 

-6.69E+03 

-8.87E+03 

-1.10E+04 

-1.32E+04 

-1.54E+04 

-1.76E+04 

-1.98E+04 



URANS, non-lineae model 
2nd= -0.01 , 4th=0.05, CFL = 1 , dt = 4^6 
0.5 M elements, 20,000 time steps 
Starting from steady RANS_NL solution 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 


LL 


n 


Vorticity Magnitude 
3.00E+04 
2.77E+04 
2.54E+04 
2.31 E+04 
2.08E+04 
1.85E+04 
1.62E+04 
1.39E+04 
1.16E+04 
9.30E+03 
7.00E+03 
4.70E+03 
2.40E+03 
1.00E+02 
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Initial condition: Linear RANS solution 


Y 


URANS, non-lineae model 
2nd= -0.01 , 4th=0.05, CFL = 1 , dt = 4e-6 
0.5 M elements, 1,5000 time steps 
Starting from steady RANS_L solution 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 




x 


I 


u 

1.32E+02 

1.21E+02 

1.09E+02 

9.71E+01 

8.54E+01 

7.37E+01 

6.20E+01 

5.03E+01 

3.85E+01 

2.68E+01 

1.51E+01 

3.42E+00 


-8.29E+00 


-2.00E+01 


Y 


URANS, non-lineae model — . 

2nd= -0.01 , 4th=0.05, CFL = 1 , dt = 4e-6 ^ x 

0.5 M elements, 1 ,5000 time steps 
Starting from steady RANS_L solution 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 


Mach 

3.01E-01 

2.80E-01 

2.58E-01 

2.37E-01 

2.15E-01 

1.94E-01 


® 1.30E-01 
1.08E-01 
8.66E-02 

■ 6.52E-02 
4.37E-02 
2.22E-02 



Y 


URANS, non-lineae model 
2nd= -0.01 , 4th=0.05, CFL = 1, dt = 4^6 
0.5 M elements, 1,5000 time steps 
Starting from steady RANS_L solution 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 



l 


turb_k 
3.79E+02 
3.52E+02 
3.25E+02 
2.98E+02 
2.71E+02 
2.44E+02 
2.17E402 
1 .90E+02 
1 .63E+02 
1 .35E+02 
1 .08E+02 
8.1 3E+01 
5.42E+01 
2.71E+01 



Y 


URANS, non-lineae model 
2nd= -0.01 , 4th=0.05, CFL = 1, dt = 4e-6 
0.5 M elements, 1,5000 time steps 
Starting from steady RANS_L solution 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 



n 


mu 

6.50E-02 
6.04E-02 
5.57E-02 
5.1 IE-02 
4.65E-02 
4.18E-02 
3.72E-02 
3.25E-02 
2.79E-02 
2.32E-02 
1.86E-02 
1.40E-02 
9.32E-03 
4.67E-03 



Y 


URANS, non-lineae model 
2nd= -0.01 , 4th=0.05, CFL = 1 , dt = 4e-6 
0.5 M elements, 1,5000 time steps 
Starting from steady RANS_L solution 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 



n 


pg 

9.30E+03 

6.84E+03 

4.37E+03 

1.91E+03 

-5.52E+02 

-3.01E+03 

-5.48E+03 

-7.94E+03 

-1.04E+04 

-1.29E+04 

-1.53E+04 

-1.78E+04 

-2.03E+04 

-2.27E+04 



ticity Magnitude 
3.00E+04 
2.77E+04 
2.54E+04 
2.31E+04 
2.08E+04 
1.85E+04 
1.62E+04 
1.39E+04 
1.16E+04 
9.30E+03 
7.00E+03 
4.70E+03 
2.40E+03 
1.00E+02 



URANS, non-lineae model 
2nd= -0.01 , 4th=0.05, CFL = 1 , dt = 4^6 
0.5 M elements, 1 ,5000 time steps 
Starting from steady RANS_L solution 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 


CL 
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Initial condition: Static flow field 


Y 


URANS, non-lineae model 
2nd= -0.01 , 4th=0.05, CFL = 1, dt = 4e-6 
0.5 M elements, 20,000 time steps 
Starting from static flow field 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 




x 


n 


u 

9.92E+01 

9.00E+01 

8.09E+01 

7.17E+01 

6.25E+01 

5.34E+01 

4.42E+01 

3.50E+01 

2.58E+01 

1.67E+01 

7.51E+00 

-1.66E+00 


-1.08E+01 


-2.00E+01 


Y 


URANS, non-lineae model 
2nd= -0.01 , 4th=0.05, CFL = 1 , dt = 4e-6 
0.5 M elements, 20,000 time steps 
Starting from static flow field 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 



I 


Mach 

3.03E-01 

2.82E-01 

2.60E-01 

2.38E-01 

2.17E-01 

1.95E-01 

1.73E-01 

1.52E-01 

1.30E-01 

1.08E-01 

8.68E-02 

6.51E-02 

4.34E-02 

2.18E-02 



Y 


URANS, non-lineae model 
2nd= -0.01 , 4th=0.05, CFL = 1 , dt = 4e-6 
0.5 M elements, 20,000 time steps 
Starting from static flow field 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 





turb_k 
1.92E402 
1 .78E+02 
1.64E402 
1 .51E+02 
1 .37E+02 
1 .23E+02 
1.10E402 
9.58E+01 
8.21E+01 
6.85E+01 
5.48E+01 
4.1 1E+01 
2.74E+01 
1.37E4C1 



Y 


URANS, non-lineae model 
2nd= -0.01 , 4th=0.05, CFL = 1 , dt = 4e-6 
0.5 M elements, 20,000 time steps 
Starting from static flow field 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 





mu 

7.42E-02 

6.89E-02 

6.36E-02 

5.83E-02 

5.30E-02 

4.77E-02 

4.24E-02 

3.71E-02 

3.18E-02 

2.65E-02 

2.12E-02 

1.59E-02 

1.06E-02 

5.33E-03 



Y 


URANS, non-lineae model 
2nd= -0.01 , 4th=0.05, CFL = 1 , dt = 4e-6 
0.5 M elements, 20,000 time steps 
Starting from static flow field 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 



I 


u 


pg 

6.76E+04 
6.52E+04 
6.28E+04 
6.04E+04 
5.80E+04 
5.56E+04 
5.33E+04 
5.09E+04 
4.85E+04 
4.61E+04 
4.37E+04 
4.1 3E+04 
3.89E+04 
3.66E+04 



URANS, non-lineae model 
2nd= -0.01 , 4th=0.05, CFL = 1 , dt = 4e-6 
0.5 M elements, 20,000 time steps 
Starting from static flow field 
Refined convective BC at exit for pg 
Extrapolated BC at exit for others 


Vorticity Magnitude 

I 3.65E+04 
3.39E+04 
3.13E+04 

2.87E+04 

2.61E+04 
2.35E+04 
2.09E+04 
1.83E+04 
1.57E+04 
1.31E+04 
1.05E+04 
7.87E+03 
5.27E+03 
2.66E+03 



NASA/TM— 2009-2 1 5486 


29 


4.1. 1.3 Centerline Variables 


The instantaneous centerline flow variables (axial velocity u , Mach number, subscale turbulent 
kinetic energy k , effective viscosity \i T , gauge pressure and vorticity magnitude) at time step 20,000 are 
presented here with three different initial conditions. From which we may examine the flow variations 
along the centerline. 
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4.1. 1.4 Convergence History 

In the present case, the time-accurate URANS simulations have reached an asymptotic steady state. 
Using the simulation in Section 4.1.1 as an example, its convergence history of the residual is given 
below. As it can be seen, the convergence of the inner (i.e., pseudo time) iteration within a physical time 
step has reached a steady (i.e., periodic) state. 


Convergence history 
URANS nonlinear model 
Unsteady convective outlet BC 



200000 400000 600000 

Iteration number 



500 r 


450 L 


400 - 


350 - 


_300 - 
PS 

2250 - 

'w 

G) 

Q£200 r 


150 - 


100 - 


Convergence history 
URANS nonlinear model 
Unsteady convective outlet BC 



774700 



J VJ 





J I I I I 1 I I I L 

774800 774900 

Iteration number 


I i 

775000 


NASA/TM — 2009-2 15486 


31 


4.1.2 Results Using Fixed Pressure at the Outlet 

Simulations have been performed using a fixed gauge pressure at the outlet while the rest of the 
dependent variables are provided by extrapolating from the interior point. The nonlinear RANS solution 
is used as the initial condition. 

The calculations have proceeded over 20,000 time steps (about 40 flow-through time) without 
difficulty. Apparently, there are very little pressure disturbances existing in the URANS solution. 

The results are presented in two parts: the time history of velocity components and gauge pressure at 
four locations along the centerline; and the instantaneous contour plots of flow variables at a center plane. 

4.1.2. 1 Time History 

The time history of velocity components and gauge pressure are recorded at four centerline locations: 
x = 0.015, 0.05, 0.10 and 0.2. From which we may examine the temporal development of the filtered 
variables. The results are presented only for the initial condition using the nonlinear RANS solution. 
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4.1.2.2 Contour at Center Plane 


The snapshots of instantaneous contour of flow variables at a center plane are presented for time step 
20,000, which is around 40 flow-through time. From which we may examine the flow structures. The 
results are presented only for the initial condition using the nonlinear RANS solution. 
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4.1.3 Results Using Extrapolation BC 

In this case, all dependent flow variables at the outlet are determined by extrapolating from the 
interior point. 

The results are presented in two parts: the time history of velocity components and gauge pressure 
along four locations of the centerline; and the instantaneous contour plots of flow variables at a center 
plane. 

4.1.3. 1 Time History 

The time histories of velocity components and gauge pressure are recorded at four centerline 
locations: x = 0.015, 0.05, 0.10 and 0.2. From which we may examine the temporal development of the 
filtered variables. The results are presented only for the initial condition using the nonlinear RANS 
solution. 
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4.1.3.2 Contour at Center Plane 

The snapshots of instantaneous contour of flow variables at a center plane are presented for time step 
of 20,000, which is about 40 flow-through time. From which we may examine the flow structures. The 
results are presented only for the initial condition using the nonlinear RANS solution. 
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4.1.4 Effect of the Outlet BC on Centerline Variables 


The above results indicate that the extrapolation boundary condition is not sustaining the dynamically 
important flow structures in the URANS simulation. This is further illustrated by inspecting the solutions 
along the centerline (see below). It can be seen that results from the unsteady convective BC and the fixed 
gauge pressure BC are almost identical; however, the results from the extrapolation BC are significantly 
different. 

Furthermore, the following results show that URANS with the unsteady convective BC (or fixed 
gauge pressure BC) leads to a sustained steady mean flow solution; however, URANS with the 
extrapolation BC can not provide a sustained mean flow solution, it keeps evolving towards a physically 
unreasonable state. 
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4.2 Linear Model 

The URANS simulations with the nonlinear subscale model are quite successful using both the 
unsteady convective outlet boundary condition and the fixed outlet gauge pressure. We have also 
performed the same simulation with a linear model (the standard k-z model) using the unsteady 
convective outlet BC. The initial condition is the nonlinear RANS solution. 

4.2.1 Results Using Unsteady Convective BC for Pressure Only 

The results are presented in two parts: the instantaneous contour plots of flow variables at a center 
plane; and the instantaneous flow variable profiles along the centerline. 

4.2.1. 1 Contour at Center Plane 

The snapshots of instantaneous contour of flow variables at a center plane are presented for time step 
of 20,000, which is about 40 flow-through time. The mean flow structures are quite similar to their 
counterparts due to the nonlinear model. 
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4.3 Comparison of Centerline Variables Between Linear and Nonlinear Models 

The flowing figures present the comparison of the distributions of the flow variables along the 
centerline between those from the linear model and those from the nonlinear model under the same initial 
condition (nonlinear RANS solution) and the same outlet boundary condition (unsteady convective BC). 
The results indicate that, there are differences in the level of the mean turbulent kinetic energy and in the 
mean vorticity distribution. 
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4.4 Concluding Remarks 

The unsteady pressure convective BC works well for URANS simulations even when the filtered 
flow is approaching to the steady state. The fixed pressure BC also works for the time accurate simulation 
of the LM6000 single injector flame tube, and its result is almost identical to that using the unsteady 
convective BC. This is not surprising, because, in this particular case, the filtered flow field itself is pretty 
steady; and there is no significant pressure disturbance near the outlet. The extrapolation BC is not 
appropriate due to its inability to sustain a physically meaningful solution. 

The initial condition has a strong effect on the evolution path of the URANS simulation towards the 
eventual state. Based on our experiences, we recommend the use of the solution from the nonlinear RANS 
simulation as the initial condition for the URANS as well as the PRNS/VLES simulations. 

Although the numerical results suggest that the URANS of linear and nonlinear k - 8 model seems to 
predict quite similar mean flow structures for the given initial and boundary conditions, there in deed 
exist appreciable differences in quantities (such as the vorticity and the turbulent kinetic energy) which 
are more sensitive to the effects of the rotation and anisotropy. 

5.0 RANS of LM6000 Single Injector Flame Tube 

The basic equations for Steady RANS approach is the Reynolds averaged Navier-Stokes equations. It 
can relatively quickly provide a global picture of the turbulent mean flow field. If the physically existing 
turbulent flow is statistically stationary, then the picture may be accurate. It can be used as a starting point 
for other higher level numerical simulations, for example, the large or the very large eddy simulation to 
explicitly bring out the unsteady turbulent structures. 

Steady RANS simulations using both nonlinear and linear models have been carried out for fixed 
pressure condition at the outlet. The purpose of these simulations is to provide the initial flow fields for 
the very large eddy simulation (PRNS/VLES) and the unsteady RANS simulation (URANS). 

5.1 Nonlinear Model 

RANS with the nonlinear model is fundamentally different from the standard k-z model. Here, the 
interactions between the turbulent mean flow and the entire spectrum of the turbulent fluctuations are 
accounted for not just by the eddy viscosity but also the explicit turbulence source terms originated from 
the nonlinear part of the model that considers the effects of anisotropy and rotation, which are missed by 
the standard k-z model. 

The results presented below are the contour plots of the mean flow variables at a center plane. It is a 
converged (over 5 orders-of-magnitude of the residual) steady solution. 

5.1.1 Contour at Center Plane 

The snapshots of mean flow variables at a center plane are presented here. It is a converged solution 
after 81374 iterations. 
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5.2 Linear Model 


RANS with the linear model is the standard k-e model. Here, the interactions between the turbulent 
mean flow and the entire spectrum of the turbulent fluctuations are accounted for by the eddy viscosity. 
Therefore, this type of model does not account for the effects of anisotropy and rotation. 

The results presented below are the contour plots of mean flow variables at a center plane. 

5.2.1 Contour at Center Plane 

The snapshots of mean flow variables at a center plane are presented here. It is a steady solution 
plotted at the 620,000 iteration. 


Y 



RANS, linear model, 2nd= -0.01 , 4th=0.05 
0.5 M elements, 620,000 iterations 


u 

n 1.30E+02 
1.19E+02 
1.08E+02 
9.74E+01 
8.67E+01 
7.60E+01 
6.52E+01 

5.45E+01 

4.38E+01 
3.31E+01 
2.24E+01 
1.17E+01 
9.53E-01 
-9.76E+00 



n_ 

RANS, linear model, 2nd= -0.01, 4th=0.05 
0.5 M elements, 620,000 iterations 




Mach 
3.10E-01 
2.88E-01 
2.66E-01 
2.44E-01 
2.22E-01 
2.00E-01 
1.78E-01 
1.55E-01 
1.33E-01 
1.1 IE-01 
8.91E-02 
6.70E-02 
4.49E-02 
2.28E-02 



Y 



RANS, linear model, 2nd= -0.01 , 4th=0.05 
0.5 M elements, 620,000 iterations 




turb_k 
4.97E+02 
4.62E+02 
4.26E+02 
3.91E+02 
3.55E+02 
3.20E+02 
2.84E+02 
2.49E+02 
2.1 3E+02 
1 .78E+02 
1.42E+02 
1 .07E+02 
7.1 0E+01 
3.55E+01 


Y 



RANS, linear model, 2nd= -0.01, 4th=0.05 
0.5 M elements, 620,000 iterations 


mu 

n 6.99E-02 
6.49E-02 
6.00E-02 
5.50E-02 
5.00E-02 
4.50E-02 
4.00E-02 
" 3.50E-02 

3.00E-02 
2.50E-02 
2.00E-02 

■ 1.50E-02 
1.00E-02 
5.02E-03 
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Y 



RANS, linear model, 2nd= -0.01, 4th=0.05 
0.5 M elements, 620,000 iterations 




pg 

6.43E+03 
4.69E+03 
2.96E+03 
1.22E+03 
-5.1 0E+02 
-2.24E+03 
-3.98E+03 
-5.71E+03 
-7.45E+03 
■9.18E+03 
-1.09E+04 
-1.26E+04 
■1.44E+04 
•1.61E+04 



£L 

RANS, linear model, 2nd= -0.01, 4th=0.05 
0.5 M elements, 620,000 iterations 


Vorticity Magnitude 

I 4.60E+04 
4.27E+04 
3.95E+04 
3.62E+04 
3.29E+04 
2.96E+04 
2.63E+04 
2.30E+04 
1.98E+04 
1.65E+04 
1.32E+04 
9.91E+03 
6.63E+03 
3.34E+03 



5.3 Comparison of Centerline Variables Between Nonlinear and Linear Models 

The following figures present the comparison of the solutions along the centerline between those 
obtained from the linear model and those from the nonlinear model. Both the contour plots and the 
centerline profiles show that the results from using the linear and the nonlinear models are quite different, 
especially near the inlet region, where the flow separation and swirling are the strongest. 
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6.0 Conclusions 

We have performed PRNS/VLES and URANS simulations for a single injector flame tube of the 
LM6000 gas turbine combustor using the National Combustion Code (NCC). Two subscale models in the 
PRNS/VLES approach and two Reynolds stress models in the URANS approach have been applied. In 
addition, in the context of each approach, three different outlet boundary conditions and three different 
initial conditions have been assessed. The major goal of these extensive numerical studies is to identify a 
relatively simple and robust outlet boundary condition for simulation approaches aimed at explicitly 
bringing out the unsteady large scale structures typically occurring in the turbulent flows of the 
combustors. 

Based on these investigations, we recommend the following outlet boundary condition: unsteady 
convective Equation (1) for the gauge pressure together with the extrapolation Equation (4) for the rest of 
the dependent variables. 

In addition, we believe that the best way to create the initial flow field for PRNS/VLES and URANS 
simulations is to run a steady RANS simulation using the nonlinear Reynolds stress model. 
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The favorite model for PRNS/VLES and URANS simulations is the nonlinear subscale model and the 

nonlinear Reynolds stress model, respectively, as the nonlinear models can account for the effects of 

rotation and anisotropy in the swirling flows. 
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